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Abstract 

The finite element method (FEM) is used to study the influence of porosity and pore shape on 
the elastic properties of model porous ceramics. The Young's modulus of each model was found to 
be practically independent of the solid Poisson's ratio. At a sufficiently high porosity, the Poisson's 
ratio of the porous models converged to a fixed value independent of the solid Poisson's ratio. The 
Young's modulus of the models is in good agreement with experimental data. We provide simple 
formulae which can be used to predict the elastic properties of ceramics, and allow the accurate 
interpretation of empirical property-porosity relations in terms of pore shape and sttucture. 

1 Introduction 

The elastic properties of two-phase (solid-pore) porous materials depend on the geometrical nature of 
the pore space and solid phase, as well as the value of porosity [|], |2|, ^ ||]. Relevant aspects of porous 
materials may include pore shape and size as well as the size and type of the interconnections between 
solid regions. These features, which generally lack precise definition, comprise the microstructure 
of the material. In order to predict properties, or properly interpret experimental property-porosity 
relationships, it is necessary to have an accurate method of relating elastic properties to porosity and 
microstructure. In this paper we use the finite element method to derive simple formulae that relate 
Young's modulus and Poisson's ratio to porosity and microstructure, for three different models of 
microstructure. 

There have been several different approaches to deriving property-porosity relations for porous 
materials. Formulae derived using the micro-mechanics method ^ ^ 0] are essentially various meth- 
ods of approximately extending exact results for small fractions of spherical or ellipsoidal pores to 
higher porosities. This includes the differential [§] and self consistent methods |l0, 11, 12] as well 



as the commonly used semi-empirical correction to the dilute result made by Coble and Kingery [[T]] to 
explain the properties of porous alumina. A drawback of this approach is that the microstructure cor- 
responding to a particular formula is not precisely known; hence agreement or disagreement with data 
can neither confirm nor reject a particular model. A second problem is that these types of models pro- 
vide no predictions for the case where the microstructure is comprised of incompletely sintered grains, 
which is a common morphology in porous ceramics. A second class of results [|3[ [l3|] have been termed 
minimum solid area (MSA) models. In this approach purely geometrical reasoning is used to predict 
the elastic moduli based on the weakest points within the structure. Again, the microstructure that 
corresponds to the MSA predictions is not exactly known. A number of semi-empirical relations have 
also been proposed [||], which generally provide a reasonable means of describing data, extrapolating 
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results and comparing data among materials. However, lacking a rigorous connection with microstruc- 



ture, these results do not offer either predictive or interpretive power. Theoretical bounds [Q, [Bp exist 
for the elastic properties, but the vanishing of the lower bound for porous materials lessens their pre- 
dictive power when the upper bound does not provide a good estimate. There are numerous other 
approaches, including the generalized method of cells [[| |J]. 

Another approach is to computationally solve the equations of elasticity for digital models of 



microstructure [J15|, |16p. In principle this can be done exactly. However, large statistical variations 
and insufficient resolution, have limited the accuracy of results obtained to date. Only recently have 
computers been able to handle the large three-dimensional models and number of computations needed 
to obtain reasonable results. As input to the method, we employ three different microstructural models 
that broadly cover the types of morphology observed in porous ceramics. The models are based 



on randomly placed spherical pores, solid spheres, and ellipsoidal pores Q14[]. The centers of the 
pores or solid particles are un-correlated which leads to realistic microstructures in which both the 
pore and solid phase are interconnected. The results, which can be expressed simply by two (or 
sometimes three) parameter relations, correspond to a particular microstructure and explicitly show 
how the properties depend on the nature of the porosity. Therefore, the results can be used as a 
predictive tool for cases where the microstructure of the ceramic is similar to one of the models, and as 
an interpretive tool if the microstructure is unknown. The numerically exact FEM results are compared 
with various well-known micro-mechanics and MSA results to determine how close an approximation 
a particular formula provides for each model. In the FEM, we can freely vary the properties of the solid 
phase, allowing us to determine the dependence of Young's modulus and Poisson's ratio on the solid 
Poisson's ratio as well as on the porosity. This question has attracted recent interest in the ceramics 
literature [jn], ||, g|j|. 



II. Computational Results 

A microstructure made up of a digital image is already naturally discretized and so lends itself to 
numerical computation of many quantities. The finite element method uses a variational formulation of 
the linear elastic equations, and finds the solution by minimizing the elastic energy via a fast conjugate 
gradient method. The digital image is assumed to have periodic boundary conditions. Details of the 
theory and programs used are reported in the papers of Garboczi & Day Jl5| ] and Garboczi [20]. 



In order to obtain accurate results using the FEM on models of random porous materials, it is 
absolutely necessary to estimate and minimize three sources of error: finite size effects, discretization 
errors, and statistical fluctuations. This has generally not been done in the past, owing to limitations 
in computer memory and speed. FEM results for random microstructures do not have much meaning 
without such an error analysis. 

The various sources of error are defined in the following way. First, the length scale of the mi- 
crostructure is fixed, usually by fixing the size of a typical pore (e.g. the spherical pore radius). The 
size of the system is then controlled by the side length of the cubic sample, denoted T. The size of 
T compared to the pore size controls how many pores will appear in the computational cell. A real 
material has many thousands or more such pores. Errors can occur in using a smaller number in a 
periodic cell to simulate a much larger number. We vary T in order to map out this effect. 

Once a value of T is chosen that minimizes finite size errors but is still computationally possible, 
we next must consider the discretization error, which comes about because we are using discrete pixels 
to represent continuum objects. The number of pixels along each edge of the cubical unit cell is M, 
giving a resolution of dx=T/M (in units of /im per pixel, if T is in /im). For the chosen value of T, 
a value of M is chosen that also gives acceptable discretization errors, usually on the order of a few 
percent. 

Finally, when computing the properties of random materials, either computationally or experimen- 
tally, one must carefully choose the number of samples (N s ) over which the results need to be averaged 
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to produce acceptable uncertainties. This value is again chosen, within computational constraints, to 
keep statistical fluctuations within a few percent. 

Overlapping Solid Spheres. Realizations of the overlapping solid sphere model [14, 21] are gener- 
ated by placing solid spheres at random points in the unit cell. This produces a set of overlapping grains 
that mimic the microstructure of sintered ceramic composites (see Fig. |]a). The space outside the solid 
grains is the pore space, with porosity 0. The pore phase is macroscopically connected above porosi- 
ties of ~ 0.03 and the solid phase remains connected for values of below 0.70 [14]. Above 
0=0.7, the solid phase is composed of isolated solid particles. So between = 0.03 and = 0.70, the 
overlapping solid sphere model is bi-continuous. In ceramics the porosity is generally less then 0.40, 
in this bi-continuous regime. We therefore consider the elastic properties for 0.1 < < 0.50, where 
the solid Poisson's ratio, u s , varied over the range 0.1 < v s < 0.4. 

To generate the microstructure we chose solid spheres of radius r=l fim. Note that the elastic 
properties are length scale invariant so the results apply to spheres of any radius for which the con- 
tinuum assumption holds. A preliminary study showed that finite size errors were acceptably small 
for cubic samples with edge length T=12 /u,m. To study the discretization errors we generated one 
realization of the model with porosity 0=0.5 at M= 48. . .128. The elastic properties depend quite 
strongly on resolution. We found that the variation of Young's modulus with M could be described 
by the relation [22] 

£fem(M) « Eo + aM- 1 (1) 

where Eq can be identified as the continuum value (corresponding to infinitely large M). The same 
is true for Poisson's ratio. Even at M = 128 the finite element code overestimates the 'exact' result 
for the Young's modulus by 30 %. Therefore, for the overlapping sphere model it is necessary to 
measure the elastic moduli at three different values of M and extrapolate the results to M — > oo. We 
chose N s =5 samples at each resolution and porosity, except at 0=0.5 where large statistical variations 
implied a larger number of samples was necessary (N S =10). Thus 30 different realizations of the 
models were considered, each at 3 different discretizations, for a total of 90 models. 

The statistical variation in Young's modulus and Poisson's ratio for the case u s =0.2 are shown 
in Table I The error bars shown in the table are equal to twice the standard error {S.E.=a/y/Ng 
with a the standard deviation). Therefore there is a 95% chance that the "true" result lies between 
the indicated error bars. The results are accurate to within 20% at 0=0.5; the error decreasing with 
porosity to less than 10% for < 0.30. The expected Gaussian distribution of the measured averages 
implies that the results are actually more accurate than this. For example, the anticipated relative errors 
are halved if a 68% likelihood threshold is used (i.e., ± one standard error). 

In addition to the above results we also computed the elastic moduli of the 90 model microstruc- 
tures at solid Poisson's ratios u s =0A, 0.3 and 0.4. The statistical variation was not significantly dif- 
ferent from the case u s =0.2. Combined with the data for v s =0.2 this covers most commonly occuring 




Figure 1: Showing pieces of the various models studied: (a) overlapping solid spheres, (b) overlap- 
ping spherical pores, and (c) overlapping ellipsoidal pores. 
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Table 1: Elastic properties of the three models (i/ s =0.2). 





Overlapping solid spheres 


Overlapping 


spherical pores 


Overlapping 


ellipsoidal pores 




E/E s 


V 




V 


E/E s 


V 


0.1 


0.71 ± 1% 


0.19 ± 1% 


0.80 ± 1% 


0.20 ± 1% 


0.73 ± 2% 


0.19 ± 3% 


0.2 


0.47 ± 2% 


0.18 ±4% 


0.62 ± 2% 


0.20 ± 2% 


0.52 ± 3% 


0.18 ±4% 


0.3 


0.25 ± 6% 


0.17 ±9% 


0.46 ±3% 


0.21 ± 3% 


0.34 ± 4% 


0.18 ± 6% 


0.4 


0.12 ± 13% 


0.15 ±25% 


0.33 ± 4% 


0.21 ± 4% 


0.20 ± 3% 


0.18 ±4% 


0.5 


0.039 ± 22% 


0.15 ±21% 


0.21 ± 8% 


0.22 ± 9% 


0.11 ±4% 


0.18 ±6% 



solids. The scaled Young's modulus for each value of u s is plotted against porosity in Fig. ||. Re- 
markably, the scaled Young's modulus of the porous material appears to be practically independent 



of u s . This result has been proven to be exact in 2-D [23, 24] and appears to hold to a very good 
approximation in 3-D. We found that the Young's modulus data are well described by an equation of 
the form 

E 



e, \ l -7»J < 2 > 

with n=2.23 and 00=0.652 and < < 0.5. Note that n and 0o are empirical correlation parameters 
and should not be interpreted as the percolation exponent and threshold, respectively. Percolation 
concepts are generally valid closer to the threshold C « 0.7 (for this model) and a higher value of n 
is expected. The computational cost of accurately measuring the elastic properties increases greatly as 
the percolation threshold is approached. 

The Poisson's ratio of the porous material is shown in Fig. |5] as a function of <j> and v s , and 



appears to be a flow diagram [23], where the Poisson's ratio asymptotically approaches a fixed point, 



independently of the value of the solid Poisson's ratio. This flow diagram has been analytically proven 



to hold in 2-D, when a percolation threshold exists at which the Young's modulus goes to zero p3|, 24]. 
This flow diagram also appears to be valid in 3-D as well, within numerical uncertainty. The Poisson's 
ratio data shown in Fig. |3| can be roughly described by the simple linear relation, 



v = v s + - (i/ - v 8 ) = Vq + 1 - - \{y s - u Q ) (3) 
<Po V <Po/ 

with two fitting parameters vq=0AA0 and O =O.472. A more accurate fit is obtained with the three 
parameter relation, 

i/ = i/ +(l-|J (v.-vq). (4) 
with v o=0. 140, 0o=O.5OO and m=1.22. 

Overlapping Spherical Pores. The overlapping spherical pore (or swiss cheese) model [J14], |2l|] is 

generated by interchanging the roles of the solid and pore phase of the overlapping solid sphere model 
(see Fig. |l|b). The morphology corresponds to isolated spherical pores at low porosity, with the pores 
becoming macroscopically interconnected at (f> 0.3. The solid phase remains connected up to « 
0.97. This type of morphology may arise in ceramics generated with a particulate filler [[I]], or where 



bubbles form in a molten state [g5J. We consider solid Poisson's ratios in the range —0.1 < v s < 0.4. 

We determined that statistical errors were acceptable for a computational cube of size T=12 /xm 
with pores of radii r=l fim. Using M=80 pixels, the discretization errors were less than 3 % for 
4> = 0.5 and 2 % for = 0.3. Therefore, it was not considered necessary to generate samples at 
different discretizations (M) and extrapolate the results. As for solid spheres, the Young's modulus 
was independent of the solid Poisson's ratio to a very good approximation. The Young's modulus can 
be described by Eq. @ with n = 1.65 and 0o=O.818 (Fig. §). Poisson's ratio of the porous material 
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Figure 2: The Young's modulus of the three microstructure models. The solid lines are empirical fits 
to the equation E/E s = (1 — cj) / 4>o) n ■ Data is shown for overlapping solid spheres (o, n = 2.23, 4>o = 
0.652), spherical pores (D,n = 1.65, <p = 0.818) and ellipsoidal pores (A,n = 2.25, <p = 0.798) 
for u s = —0.1, . . . , 0.4. Note that the E is practically independent of the solid Poisson's ratio in each 
case [the different values of E[y s ) at each porosity are almost indistinguishable]. 




Figure 3: The Poisson's ratio of the overlapping solid sphere model as a function of porosity for 
v s =0. 1-0.4. The dashed lines are an empirical fit to Eq. (|3j> . The solid lines correspond to the three 
parameter relation given in Eq. (||), with the value of all parameters given in the text. The intercepts 
of the lines at zero porosity correspond to the solid Poisson's ratio. 



is shown in Fig. |] and is simply described by the linear relation given in Eq. (||) with v$=0.22\ and 
^0=0-840. Again, a flow diagram is observed. 

Overlapping Ellipsoidal Pores. A common method of analyzing the effect of pore shape on elastic 
properties is to study ellipsoidal pores. In analytic formulae, it is possible to treat the limiting cases 
of needles and platelets, although the difficulty of resolving these fine structures prohibits these lim- 
its from being treated with the finite element method. However, the percolation properties of these 
limiting cases can be computationally studied [26]. To gauge the effect of deviations from spherical 
shaped pores we considered isotropically oriented overlapping oblate ellipsoidal pores bounded by 
the surface (x/a) 2 + (y/b) 2 + (z/c) 2 = 1 with a=b=l fim and c=0.25 ^m (see Fig. [j]c). For this 
case, the pore phase becomes connected at porosity <ft = 0.2 [|2(J]. Statistical errors were found to be 
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Figure 4: Poisson's ratio of the overlapping spherical pore model as a function of solid Poisson's ratio 
and porosity. The lines are an empirical fit to the relation v = v s + <p/<po x (z/ — v s ) with ^=0.22 1 
and (/>o=0.840. The intercepts of the lines at zero porosity correspond to the solid Poisson's ratio. 




0.0 0.2 0.4 0.6 

Porosity 

Figure 5: Poisson's ratio of the overlapping ellipsoidal pore model as a function of solid Poisson's 
ratio and porosity. The solid lines are an empirical fit to Eq. (Q) and the dashed lines correspond to 
the linear fit to Eq. (|3|) with uq = 0.166 and 4>q = 0.604. The intercepts of the lines at zero porosity 
correspond to the solid Poisson's ratio. 

acceptable for a computational cube of size T=10 /mi. Using M=96 pixels the discretization errors 
were 3 % for <fi = 0.5 and 2 % for (ft = 0.3. As for the case of spherical pores, these errors were 
considered sufficiently small, so that the added computational burden of the extrapolation technique 
could be again avoided. 

Again the Young's modulus was found to be independent of the solid Poisson's ratio to a very 
good approximation. The results, shown in Fig. ||, can be accurately described by Eq. @ with n=2.25 
and 0o = O-798. The Poisson's ratio (Fig. |5|) can be roughly fit using the Eq. (||) with i^)=0.166 and 
00=0.604. A better fit is obtained using Eq. @ with m=\.9\, z/ =0.161 and O =O.959. A flow diagram 
similar to that seen before is obtained. 

The CPU time and memory required for these computations are an important "experimental" de- 
tail. The memory requirement for a given model was 230 x M 3 bytes, where M was the edge length 
in pixels of a cubic unit cell. So for the largest computations carried out, M = 128, the memory 
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requirement was about 500 Mbytes. The amount of CPU time consumed was approximately 3000 
hours, divided among different modern workstations. 



III. Comparison with micro-mechanical and MSA formulae 

In this section we compare a selection of well known theoretical results with the 'numerically exact' 
data computed in the previous section. These results include both analytically exact results (bounds, 
expansions, dilute limits, composite sphere assemblage), as well as approximate results (effective 
medium theories, minimum solid area models). 

There are several kinds of exact bounds that have been derived for elastic materials [Q, [I4j]. These 
are equations involving the different phase moduli, the volume fractions of the different phases, and 
various correlation functions that define the geometry of the composite. The upper bound gives the 
maximum possible composite elastic moduli, and the lower bound gives the lowest possible composite 
elastic moduli. The bounds used in this paper are three- point bounds, which have been written out 
explicitly for overlapping solid spheres and overlapping spherical pores [Q]. In the case where one 
phase has zero elastic moduli, as is true in this paper, the lower bound becomes zero as well, and so 
only the upper bound is meaningful. 

An exact perturbation expansion also exists, where the elastic moduli of a two-phase material is 
expanded in terms of parameters involving the individual elastic moduli of each phase and geometrical 



quantities [ g7j [28]]. This expansion has been carried out to three terms explicitly, and it is this truncated 
form to which we will compare our numerical data. The result is expected to be accurate when the 
void phase is not interconnected. 

Another exact result, which is used later in this section to build the various effective medium 
theories, is the case of dilute spherical pores for which the exact effective moduli are given by 

K = K m + c i P mi {K i -K m ) (5) 
G = G m + Cl Q mi (Gi-G m ) (6) 



with 



3K m • \G at (^"" Gin H ■ 1 1, 



3Ki + 4G m ^ d + F, 
G m 9K m + 8G 
6 K m -\- 2>G m 



(7) 



Here Cj denotes the concentration (volume fraction) of inclusions and the subscripts i and m on the 
bulk K and shear modulus G denote the properties of the inclusion and matrix, respectively. The result 
is attributed to numerous authors [Q]. For a porous matrix Ki = G{ = and the porosity is (j) = a. 
The result is strictly valid for small concentrations of inclusions <p <C 1 (in practice < 0.1.) Cast in 
terms of the engineering constants for porous inclusions this result becomes 

E = Em -hE j~^ m ~ 5 ^ +O(0 2 ) (8) 
2 7 - bum 

„ _ 3 ^.,.-ixi-,;,) 2 

2 7 - bu m 



Our prior statement p9Q of Eq. (JS|> inadvertently omitted the factor of 3/2, although the correct 
result was used in the paper. A non zero quadratic term can be added (as an empirical correction) to 
ensure that E=0 at 0=1. This was suggested by Coble and Kingery [|T|] for MacKenzie's [ [30| ] result 
for spherical pores, which is equivalent to Eqs. with K{ = G{ = 0. 

To adapt the dilute formulas to the case of finite porosity a number of proposals have been made. 
The approximate equations that result are usually called effective medium theories. The most com- 



mon approximation is the so-called self consistent method (SCM) of Hill [|9|] and Budiansky QIOQ. In 
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this model the equations of elasticity are solved for a spherical inclusion embedded in a medium of 
unknown effective moduli. The effective moduli K and G are then derived. In the dilute case the 



embedding medium is just the matrix. The Hill-Budiansky result can be stated as [ 12] 



dP* i {K i -K )f ) + c m P* rn {K m -K )f ) = 
CiQ* l (Gi — G*) + c m Q* m (G m — G^) = 



(10) 
(11) 



where K* and G* denote the effective moduli and P* m and Q* m are given in Eq. (^). The equations 
cannot be explicitly solved and numerical methods are necessary (see Hill [§] and Berryman [12] 
for details). In the case of porous inclusions, the moduli vanish at (ft = |, which is a property not 
shared with most composites (e.g., the overlapping sphere model). To derive a more realistic result, 



Christensen and Lo [ pip generalized the SCM (GSCM) to the case of a spherical shell embedded in a 
matrix of unknown moduli. The result is complicated and not reproduced here. 

The differential method (reviewed by McLaughlin [Jsj] ) provides an alternative model using a simi- 
lar philosophy. Suppose that the effective moduli of a composite medium are known to be K+ and G*. 
Now if a small additional concentration of inclusions are added, the change in and G* is approx- 
imated to be that which would arise if a dilute concentration of inclusions were added to a uniform, 
homogeneous matrix with moduli and G* . This leads to a pair of coupled differential equations, 



dK* 
da 
dG* 
da 



; Ki — K. 



P 



1 - Cj 

Gj — G 



-; K*{c i = 0)=K r> 
; G*(cj = 0) = G m 



(12) 
(13) 

] can all be extended 
have been given by 

Berryman Q12Q. In addition to these results Wu Jll[ ] derived a variant of the self consistent method, 
where K* and G*, the effective moduli, are found by implicitly solving the equations 



The dilute result, the self consistent result [12Q, and the differential method | 
to the case of spheroidal inclusions. The general results [11] for P m% and Q r 



= K m + c i P* i {K i -K m ) 
G* = G m + CiQ* l (Gi — G m ). 



(14) 
(15) 



A different type of microstructure is provided by Hashin's [|32|, ||] model of space-filling poly- 
disperse hollow spheres (the "composite-sphere assemblage"). Although a simple formula exists for 



the bulk modulus over the full porosity range [32], exact results for the Young's moduli are not avail- 



able. Ramakrishnan and Arunachalam [ |33[ ] recently derived the approximation 



E 
El 



(1 



(1 + 2<f> - 3v s <p) 
(&v s + 3(f) - 7v s c, 



4(1 + 



3za, 



(16) 
(17) 



However, the derivation is not rigorous. In particular, the exact result for the bulk modulus of the 

0.5. Since Eq. ( fl6| ) was found 



model Q32[ ] is around twice that predicted by Eqs. (|16|-|17j) at </ 

to provide reasonable agreement with experimental data for porous ceramics [|l7j|, we compare its 
predictions to our FEM data below. 



The final class of results we consider is provided by the 'minimum solid area' (MSA) models [ 34 ] 
(which have been recently reviewed by Rice |l3]]). This approach is based on the assumption that 
the ratio of the effective moduli to the solid moduli is directly proportional to the minimum ratio of 
solid contact area to the total cross-sectional area of periodic structures. The approximation derived 
depends on the particular model considered. We consider two basic models most closely aligned with 
our FEM data: simple cubic arrays of solid and porous spheres. The latter case provides a particularly 



8 




o.o H 1 1 1 

0.0 0.2 0.4 0.6 

Porosity 



Figure 6: A comparison of rigorous bounds and expansions to the FEM data for overlapping spherical 

pores (o) and overlapping solid spheres (□). The truncated expansion ( ) and the three-point upper 

bound ( ) are shown for the spherical pore case. Only the three-point bound (• • •) is shown for the 

solid sphere case. The three-point lower bound is zero for porous materials. The Poisson's ratio is 0.2 
for all the results. 



simple example of the type of result which can be derived. Suppose the repeat distance of the lattice is 
2h and the sphere radius is r. The Young's modulus is assumed to be proportional to the area fraction, 
giving 



2 



E (2h) 2 -irr 2 vr /6\f ,2 

- 1 - - - W>3 (18) 



E s (2h) 2 4 Vvr 

since cf> = ^ir(r/h) 3 . The form of the result changes for r > h (or (f> > 7r/6=0.52) as the spheres begin 
to coalesce. Rice [||] has noted the results of many different periodic structures can be approximated 
by the form E/E 3 = e~ b< ^ over a range of porosities. For example, 6 5 for the solid sphere model 
and b rs 3 for the porous sphere model. It is argued that for a given set of data, b can be compared with 
known values to assess the type of porosity. Often fractions of different types of porosity are assumed 
to match experimental data making the method an interpretive rather than a predictive tool. Since we 
have measured E for microstructures based on solid sphere contacts and porous spheres we should be 
able to ascertain the accuracy of the MSA formulae for these cases. 

Fig. [] shows the comparison between the exact three-point bounds [ 14 ] for the overlapping solid 



sphere and spherical pore case, the truncated expansion [ |27| , |28[ ] for the overlapping spherical pore 
case, and the numerical results. Clearly the expansion does better than the three-point bound for the 
overlapping spherical pore case, though both are fairly close to the numerical results. The bound lies 
far away from the overlapping solid sphere numerical results, however. For this case, the truncated 
expansion does not exist. Only the v s = 0.2 data is shown. Using the truncated expansion, one can 
show that in 3-D, the Young's modulus is not exactly independent of the solid Poisson's ratio, but is 
rather a very good approximation, as was shown earlier in this paper. 

In Fig. ^ we compare the FEM data (y s = 0.2) for overlapping spherical pores with dilute and 
effective medium theory analytic results. At this Poisson's ratio the SCM and dilute results reduce to 
E/E s = 1 — 2<f) while the differential and dilute results with the Coble-Kingery correction reduce 
to E/E s = (1 — (p) 2 . Since the analytic results are based on the case of dilute spherical pores they 
all match the FEM data at = 0.1. The dilute and SCM results under-estimate the FEM data at 
higher porosities because of the aphysical percolation threshold at <fi = |. The generalized SCM over- 
estimates the data, while the differential method performs reasonably well over the entire porosity 
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Figure 7: A comparison of different theories to the FEM data for overlapping spherical pores (o). The 

lines correspond to the dilute result and self consistent method [[J [K]] (or SCM) ( ), the differential 

method [||] and dilute result with Coble-Kingery correction ( ) and the generalized SCM [31] ( — 



). Data for the overlapping solid sphere model (for which no rigorous theories exist) are also 

shown (□). 



range. The latter observation might have been anticipated given the close association between the 
definition of the spherical pore model and the assumption of the differential method. At increasing 
porosities we are simply adding additional spherical pores to a porous matrix. The data for overlapping 
solid spheres is also shown in the figure, and seen to be quite different from any of the available results. 
This demonstrates that microstructure (the geometrical nature of the porosity) is an important factor 
besides the actual value of the porosity 

In Fig. ||, the minimum solid area models and the Ramakrishnan and Arunachalam results [33] 
are compared with the data. The MSA model for spherical pores performs reasonably well, although 
underestimating the FEM data for overlapping random spherical pores at low porosities <ft < 0.3. The 
MSA model for solid spheres considerably underestimates these data for (p < 0.3. The Ramakrish- 
nan and Arunachalam [ f33| ] approximation falls between the FEM data for <f> > 0.1 indicating that it 
corresponds to neither of the microstructures. For purposes of comparison we also report numerical 
results obtained using the computational generalized method of cells [Q]. For aperiodic spherical pore 
the results significantly underestimate the FEM data for overlapping spherical pores at low porosities 
(and hence the exact dilute result). It is not clear if this is due to the assumptions, or the particular 
implementation, of the method. 

The FEM data for overlapping oblate ellipsoidal pores is compared with the available theories 
in Fig. |9| The SCM results of Wu [11] and Berryman [ |l2| ] underestimate the porosity as a result 
of underestimating the physical percolation threshold. The Berryman result performs significantly 
better than does the Wu result. As for the case of spheres, the differential method matches the data 
quite closely because of the similarity between the assumptions of the theory and the definition of the 
model. 

We have also compared the Poisson's ratio predicted by the various self-consistent and differential 
methods to the FEM data for overlapping spherical and ellipsoidal pores. The theoretical results 
converge to different fixed points (e.g. Fig. |3]) in qualitative agreement with the data. But only the 
differential method provides reasonable agreement with the FEM data (with absolute error less than 
0.02 for (p < 0.4 and 0.1 < u s < 0.4). 



IV. Comparison with Experiment 
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Figure 8: A comparison of the minimum solid area (MSA) models [^] to the FEM data for overlapping 

spherical pores (o) and solid spheres (□). The MSA solid sphere model ( ) and MSA porous 

sphere model ( ) (in simple cubic packings) are seen to under-estimate the data for low porosities 

(4> < 0.3). The formula of Ramakrishnan and Arunachalam [Q E/E s = (1 - <j)) 2 /(l + 1-40) (• • •) 
and the results of the generalized method of cells for a periodic spherical pore [§](-•-) are also 
shown. 




Figure 9: A comparison of different theories to the FEM data for overlapping oblate ellipsoidal pores 

(o). The lines correspond to the differential method [||] ( ) and the self consistent methods of 

Wd ifTh ( ) and Berryman [O] (• • •)■ 



We now use the FEM results to analyze experimental measurements of the elastic properties of porous 
ceramic materials. The dependence of the elastic moduli on porosity has been the subject of many 
studies [13, 17|, pH| ]. Data for porous alumina from numerous studies [^] are shown in Fig. [H]. The 
Coble-Kingery ^ material is markedly stiffer than other materials, and is in very good agreement 
with the FEM results for the overlapping spherical pore model. The pores in the alumina matrix 
were actually created by the incorporation of a particulate filler [pi], which corresponds well with the 
definition of the model microstructure. The remaining data closely follow the overlapping solid sphere 
FEM result for <fi < 0.25, indicating that the solid alumina phase has the sintered granular morphology 
exhibited by the model microstructure (Figjlja). However, Knudsen notes that several of the samples 
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Figure 10: Data for alumina (E S =410 GPa) compiled by Knudsen [35] (o). The Coble-Kingery [|]] 
(£' s =386 GPa) data are also shown (□). The lines correspond to the FEM theories computed in this 

paper: overlapping spherical pores ( ), overlapping oblate ellipsoidal pores (• • •) and overlapping 

solid spheres ( ). 



summarized were also created using particulate fillers. At higher porosities the solid sphere result 
underestimates the data. One reason for this might be that the model contains isolated solid spheres 
which artificially reduce the actual porosity. This was checked and found not to be the case for the 
porosities studied. Therefore, the solid connections in these samples of porous alumina are likely 
stiffer than those found in the solid sphere model at porosities <ft > 0.25. Overlapping spheres can 
create very sharp "valleys" between a pair of overlapping solid spheres (see Fig. |I]a), which would be 
rounded off in the sintering process, presumably strengthening the solid-solid connection. 

Hunter et al. [36, [3^, 38, 39] have studied the Young's modulus of several different oxides. In all 
cases, the porous material was created by sintering a powder of the pure oxide. The results for the 
Young's modulus are reproduced in Fig. 11. For low porosities {<j) < 0.1) all of the data followed the 
FEM results for overlapping spherical pores. For Gd203 the FEM result continues to provide excellent 
agreement up to the maximum porosity measured (0 = 0.4) indicating that the microstructure is sim- 
ilar to that of the model (overlapping pores). In contrast, the data for the other three oxides decreases 
towards the result for overlapping solid spheres indicating a more granular character. 

The data of Walsh et al. [25] for porous glass is compared with the FEM results for overlapping 
spherical pores in Fig. 12. The agreement is good for small to moderate porosities (</> < 0.3), but the 
FEM results underestimate the data at higher porosities. Walsh et al. point out that the pores in the 
glass are actually not interconnected (unlike the overlapping pores of the model). This would account 
for the increased stiffness. It is interesting that the FEM results begin to deviate from the experimental 
data at the threshold where the pores become macroscopically connected (0=0.3). Data for sintered 
MgAl202 [ |40| ] powder is shown in Fig. 13, and is well modeled by the FEM results for overlapping 
solid spheres. Micrographs of the ceramic indicate a granular structure similar to that of the model 
microstructure (although the grains appear more like polyhedra, not spheres). 
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Figure 11: Data for various oxides measured by Hunter et al [3^, 37, 38, 3^] compared with the FEM 

theories for overlapping spherical pores ( ), and overlapping solid spheres ( ). Sni203 [ f36| ] 

E S =U5 GPa (o); Lu 2 3 @] E s =\93 GPa (A); Gd 2 3 @ £ s =150 GPa (o); Hf0 2 @ E 8 =2A6 GPa 

(□); 




Figure 12: Data for porous glass [ ]25p (K s =46 GPa, v s =0.23). The line corresponds to the FEM theory 
for overlapping spherical pores ( ). 
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Figure 13: Data for MgAl 2 4 [@]. We used the value E s =4l.2x 10 6 psi C 2 ^ 4 GPa) indicated on 
Fig. 3(A) of the reference, rather than the reported value of £' s =43.4xl0 6 psi which appears to be a 
misprint. The line corresponds to the FEM theory for overlapping solid spheres. 



V. Discussion and Conclusions 

We have derived empirical theories for the dependence of the Young's modulus on porosity for three 
distinct models of porous ceramics, based on careful finite element computations. An advantage of 
these results over many conventional theories is that they correspond to a priori known microstruc- 
tures. The dilute result (extended by Coble-Kingery to all porosities), the differential method, and 
the self consistent methods all have a "built-in" microstructure, but apart from the dilute case, it is 
not clear what that microstructure is. Therefore, agreement or disagreement with a particular analytic 
theory neither confirms nor rejects a particular physical model. For the minimum solid area models, 
the microstructure is exactly known, but the approximation involved in making the Young's modulus 
directly proportional to the contact area leads to a similar conclusion. Indeed, we found that the MSA 
models did not provide quantitative agreement with the moduli of the random microstructures stud- 
ied. We found that the differential method (Eqs. [l2| and |l3|) gave results in reasonable agreement with 
computed data for the cases of overlapping spherical and ellipsoidal pores, probably due to the simi- 
larities between the assumptions of the model and the definitions of the microstructure. Results for the 
granular model of overlapping solid spheres were not well modeled by any of the analytic theories, 
demonstrating the importance of finite element techniques in this case of great physical interest. 

We have also generated data that shows the dependence of Poisson's ratio on porosity and the solid 
Poisson's ratio. It is difficult to study this question experimentally because of the inability to vary the 
Poisson's ratio of the solid independently, and the well known difficulties of accurately measuring 



the Poisson's ratio at moderate to high porosities [|18p. At sufficiently high porosities we find that 
the Poisson's ratio converges to a fixed non-zero value uq irrespective of the solid Poisson's ratio. 
For overlapping solid spheres uq = 0.14, spherical pores uq = 0.22 and oblate ellipsoidal pores 



vq = 0.16. This behavior is exact in two dimensions [123 , 24] and is exhibited by many of the analytic 



theories in three dimensions. At present the available experimental data cannot confirm this qualitative 



behavior []18[]. We have shown that the Poisson's ratio does not vanish at high porosities as has been 
recently argued [|T9|]. 

It is not simple to attribute our results to features of the solid-pore morphology such as the size, 
shape, distribution and connectivity of pores or solid grains, since these features have no obvious 
definition for complex bi-continuous random microstructures. A few general observations can be 
made, and interpreted in terms of interrelated geometrical and mechanical features of the models. For 
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a given porosity, the sintered grain structure of the overlapping solid sphere model is relatively weak. 
The small solid contacts between spheres and the highly interconnected porosity (which becomes 
macroscopically connected at 0=0.03) lead to a weak structure. We also assume that the valleys 
which occur between grains will provide sites of large stress concentrations, and consequently, large 
deformations. In contrast, spherical pores provide high (near optimal) stiffness at a given porosity. The 
dispersed nature of the porosity (which is macroscopically disconnected for < 0.3) corresponds to a 
well connected solid matrix. Ellipsoidal pores tend to weaken a structure more than spherical pores due 
to a combination of a less well connected solid phase (the pores become macroscopically connected 
at <p = 0.2), and greater stresses and deformations near the high curvature regions of the ellipsoid. 

We have compared our FEM results with several sets of previously published experimental data. In 
cases where the microstructure of the porous ceramics roughly matched that of the models, the agree- 
ment was very good. Since the FEM results correspond to a known microstructure, it was possible 
to explain deviations in terms of specific microstructural features. Thus, comparison of experimental 
data with the three computational results provides a useful interpretive tool. Note that a given elastic 
modulus does not correspond to a particular microstructure. Therefore, it is important to corroborate 
microstructural interpretations obtained from the elastic moduli with information about the particu- 
lar material (such as a micrograph). In the future it would be useful to extend this work to higher 
porosities and to other relevant models (such as non-overlapping porous spheres). It is also possible 
to use statistical microstructural information obtained from two-dimensional micrographs to generate 
models p9| ] that actually mimic physical microstructures. 
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